Mode entanglement of electrons in the one-dimensional Prenkel-Kontorova model 
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We study the mode entanglement in the one-dimensional Frenkel-Kontorova model, and found 
that behaviors quantum entanglement are distinct before and after the transition by breaking of 
analyticity. We show that the more extended the electron is, the more entangled the correspond- 
ing state. Finally, a quantitative relation is given between the average square of the concurrence 
quantifying the degree of entanglement and the participation ratio characterizing the degree of 
localization. 
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Quantum information science have emerged as an ac- 
tive interdisciplinary area between quantum mechanics 
and information theory 0. Recently, it was suggested 
that the quantum information science may offer an pow- 
erful approach to the study of nonlinear complex quan- 
tum systems |^ |E Specifically, there may be close 
connections between quantum entanglement theory and 
many-body theory P, 0]. These ideas motivate us to 
investigate nonlinear complex quantum systems by en- 
tanglement theory. 

In this letter, we study quantum entanglement of 
electrons in the Frenkel-Kontorova (FK) model a 
paradigm in nonlinear science, and address the effects 
of the transition by breaking of analyticity (one striking 
feature of the FK model) on behaviors of entanglement. 
The FK model has been used to model various kinds of 
physical systems such as an electron in a quasi-lD metal 
below the Peierls transition ^] . The FK model describes 
a one-dimensional chain of atoms with harmonic nearest 
neighbor interaction placed in a periodic potential. Due 
to the competition between the two length scales, the 
spring length and the period of the on-site potential, the 
FK model exhibits a rich complex phenomena -"-^ 
shown by Aubry that there exist two different ground- 
state configurations for an incommensurate chain, and 
the transition from one configuration to another is driven 
by a single parameter K. These two incommensurate 
configurations correspond to invariance circle and can- 
torus of the standard map 0| , respectively. 

The electronic properties such as the energy spectrum 
and quantum diffusion in the FK model have been stud- 
ied |lll|. Quite recently, entanglement properties were 



studied in the Harper model |12| . another paradigm of 
nonlinear science, and some connections are revealed be- 
tween entanglement and localization [l^ . In contrast to 
the Harper model that have been often used to s tudy 
electron properties in incommensurate systems |l2L Il4j , 
the FK model has two control parameters K and A, 



which leads to more rich physics. Moreover, the Harper 
model exhibits a symmetry of self-duality, whereas the 
FK model does not. Next, we study ground-state entan- 
glement properties of the FK model, and find that the 
entanglement changes drastically when one goes from one 
configuration of atoms to another. 

Let us start by recalling some basic facts about (spin- 
less) fermions on a lattice of atoms. Consider N lo- 
cal ferniionic modes (LFMs) ~ sites which can be either 
empty or occupied by an electron [isj . In the second- 
quantized picture the basic objects are the creation and 
annihilation operators cjj and c„ of n-th LFM, satisfying 
the canonical anti-commutation relations 
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The Hilbert space naturally associated to the N LFMs, 
known as Fock space 7i f , is spanned by 2^ basis vectors 
|ni,...,ni) :=n/=i(c|)"'|0> (n, = 0,lVO. 

From the above occupation-number basis it should be 
evident that Tip isomorphic to the A^-qubit space. 
This is easily seen by defining the mapping [iq 
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A:=Y[{clr\0)^®l,\ni) 
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where a'^ is the raising operator of ^-th qubit. This is 
a Hilbert-space isomorphism between Hp and C®^. Af- 
ter this identification one can discuss entanglement of 
fermions by studying entanglement of qubits. Clearly, 
this entanglement is relative to the mapping JSJ. By 
defining new fermionic modes via automorphisms of the 
algebra (Q) one gives rise to different mappings between 
Hp and with an associated different entanglement. 
This simple fact is one of the manifestations of the rela- 
tivity of entanglement 

Consider an electron hopping in a ID FK chain de- 
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scribed by the following Hamiltonian 
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CnCn, 
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where t is a nearest-neighbor hopping integral which is set 
to 1 throughout the paper, and Vn — A cos(27r(Ta;° ) is the 
on-site potential, which is controlled by the parameter A 
and the configuration {x^}. Here, A is the amplitude of 
the on-site potential, a = Fn/Fn-i is the inverse distance 
between two consecutive atoms for K = 0, and {-Fk} is 
a Fibonacci sequence and the series of truncated fraction 
Fn-i/Fn converse to the inverse golden mean (\/5 — 1) /2. 
The number of lattice sites is chosen to be = F„. 

The configuration {x^} is the FK model is determined 
by minimizing the functional 

U = Y^ ^{xn+i - Xn f + K[l - cos(27ra;„)], (4) 

n 

where K is a coupling constant controlling configurations 
of atoms. Moreover, the periodic boundary condition is 
assumed for Hamiltonian 

A general state of a single electron moving on the ID 
FK chain is given by 
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where |0) is the vacuum state. From Eqs. ijSJl and (O, we 
obtain the eigenequation of the system 



(6) 



from which we can obtain the ground state. One-particle 
ground states are always extended for periodic systems, 
which can be localized for the FK model as the on-site 
potential Vn may lead to aperiodicity. 

To study entanglement properties, we first map a state 
of fermions to that of spins according to Eq. |2Jl. Af- 
ter the mapping, our state l^f) becomes a multi-qubit 
state, and the pairwise entanglement quantified by the 
concurrence [l^ is well-defined. For state (Q, the con- 
currence for two LFMs i and j is easily found to be 
Cij = 2|i/)iV'j|- Specifically, when Itj^nl — ^/VN, the state 
becomes the so-called W state |2y and the concurrence 
is given by 2/N. In this study, we are more interested in 
the gross measure of entanglement, the average concur- 
rence nal 
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which have connections to localization. Here, M ~ 
TV (TV — l)/2. We will mainly concentrate the pairwise 
entanglement of ground states and briefly discuss the bi- 
partite entanglement. 




FIG. 1: Average concurrence against parameter K for FK 
chains with different lengths. The parameter A = 3. 



As a first step of numerical calculations, we obtain the 
configuration for N atoms by gradient method , adopt- 
ing the periodic boundary condition, ft is well-known 
that there exists a critical value Kc = 0.154641 sepa- 
rating two configurations of atoms. The configurations 
determine the on-site potential and thus ground-state 
entanglement. 

Figure 1 displays behaviors of the average concurrence 
of the ground state of the FK model as a function 
of K. When the parameter K increases, we observe an 
abrupt decrease of the average concurrence near criti- 
cal value Kc- There is a strong interrelation between 
the electronic properties and the configurations of atoms, 
and thus this transition of entanglement results from the 
transition by breaking of analyticity in the configurations 
of atoms. The entanglement strongly feels the classical 
transition. For K < Kc, the configuration of atoms cor- 
responds to invariance circle. In this case, the concur- 
rences (C) w 2/iV, and the ground states are extended. 
For K > Kc, the configuration corresponds to cantorus, 
and ground states are quite different from the case of 
K < Kc- In this case, the concurrence tends to disap- 
pear and the ground state is localized. 

To see clearly between the entanglement and localiza- 
tion, we show in Fig. 2 the average concurrence against 
the participation ratio. The participation ratio charac- 
terizing the degree of localization is defined by 



1 



P 



(8) 



We see that the average concurrence increases with the 
increase of the participation ratio, illustrating that the 
more localized ground states are, the less the average 
pairwise entanglement. 

Now, we investigate effects of amplitude A of the on- 
site potential on ground-state entanglement. The average 
concurrence as a function of A is shown in Fig. 3 for dif- 
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FIG. 2: Average concurrence against the participation ratio. 
The parameters A = 3 and TV = 377. 
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FIG. 3: Average concurrence against the parameter A for dif- 
ferent K. The parameter A*' = 377. 



ferent K . When K = 0.1, the concurrence is nearly un- 
changed, and approximately given by 2/377 = 0.005305. 
As if < Kc, atoms are in the configuration before the 
breaking of analyticity, and this configuration is nearly 
the same as that for K = 0. Then, the on-site potential 
Vn = A cos(27r(Ta;° ) are approximately of no difference, 
and the second term of Eq. Q contribute a constant 
to the Hamiltonian for one-particle states. The concur- 
rence of the ground state is then only determined by the 
first term of Eq. (O, and thus independent of A. When 
K = 0.3, the concurrence reduces quickly to zero as A 
increases. In this case, the atoms are in the configu- 
ration after the breaking of analyticity, and the second 
term of Eq. (PJ no longer commute with the first one 
and has significant effects on the ground-state proper- 
ties. For larger A, the second term dominate over the 
first one, and thus reduce the entanglement. It is evident 
that the concurrence goes to zero as A ^ oo. The curves 
for JsT = 0.16 and 0.165 displays the intermediate behav- 
iors. For K = 0.3, we observe a critical value Ac, after 
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FIG. 4: Time evolution of the average concurrence. The pa- 
rameter N = 377. (a), (b) and (c) correspond to A = 1, A = 2, 
and A = 4, respectively; K = 0.1. (d) corresponds to K = 0.3 
for different A. The parameter N = 233. 



which the entanglement tends to zero. The critical value 
is dependent on the parameter K. 

So far, we have studied ground-state entanglement 
properties in the FK model. Next, we investigate dy- 
namics of entanglement with the initial state being |1) = 
c||0). Thus, there is no initial entanglement. The 
time evolution is described by a time-dependent equa- 
tion {h = 1) 



.dtlJn , 



+ 1 



i^n-i + Xcos{2nax1)^pn, (9) 



which can be integrated numerically by various methods 
such as the fourth-order Runge-Kutta method. 

Figure 4 displays dynamical behaviors of the average 
concurrence for different K and A. From Fig. 1 (a), (b) 
and (c) for if = 0.1, we observe that the concurrence first 
linearly increases with time, and finally reaches a plateau, 
oscillating irregularly around a steady value. The dy- 
namical behaviors only differ slightly for different A. The 
reason is similar to that discussed above, namely, for this 
case the term A cos(27r(TX° ) contributes a constant to the 
Hamiltonian and does not alter the system dynamics. In 
contrast, for K = 0.3, as seen from Fig. 1 (d), the increase 
of A suppresses the entanglement generation. When A 
becomes larger, the electronic states becomes more lo- 
calized, and the entanglement diminishes. The slight de- 
pendence of A for if = 0.1 and strong dependence of A for 
K — 0.3 corresponds to extended and locahzed ground 
states, respectively. 

From the above results on ground-state and dynami- 
cal behaviors of entanglement, we see that the more lo- 
calized a state is, the less the entanglement. Now, we 
build a direct connection between the concurrence, quan- 
tifying the pairwise entanglement, and the participation 
ratio, characterizing the degree of localization. For one- 
particle state I^E') (0, the concurrence between LFMs i 
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and j Cij = 2\tpiipj\. We make an average of square of 
concurrence, rather than the concurrence, 

i<j i<j 
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Np 



(10) 



In deriving of the above equation, we have used Eq. ^ 
and the identity 2 X;i<j I'/'iPlV'jP = ^-J2iLi IV'il'*, which 
results form the normahzation condition X)t=i ~ ^■ 
Thus, the average of the square of concurrence can be 
written as a simple function of the participation ratio, 
and this relation build a direct connection between pair- 
wise entanglement and localization. It is evident that 
the larger p is, the larger the concurrence. For the two 
extreme cases, p = 1/N and p — 1, the concurrence C 
becomes and A/N"^ , respectively, as we expected. Note 
that the relation is applicable to arbitrary one-particle 
states, irrespective of model Hamiltonians. 

We have discussed the pairwise entanglement. For 
other type of entanglement, such as the bipartite pure- 
state entanglement quantified by the linear entropy, we 
also find similar relations as Eq. (|l()ll between the lin- 
ear entropy and participation ratio 21]. The connec- 
tions between entanglement and localization are not re- 
stricted to electronic system, and can be applied to other 
systems such as spin systems with one-magnon excita- 
tions. The investigation of multipartite entanglement 
other than pairwise and bipartite entanglement is more 
interesting, but at the same time more difficult and com- 
plicated. 

In conclusion, we have studied ground-state and dy- 
namical pairwise entanglement of two LFMs in the one- 
dimensional FK model, and found that the entanglement 
exhibits distinct behaviors for the cases oi K < Kc and 
K > Kc- This is a consequence of the transition by 
breaking of analyticity. For K < Kc, the ground state 
is extended and more entangled; while for K > Kc the 
ground state is localized and less or not entangled. The 
amplitude A of the on-site potential have slight effects on 
entanglement when K < Kc, while has significant effects 
when K > Kc- It is interesting to note that entangle- 
ment is closely connected to localization. It becomes a 
general feature that the more extended the electron is, 
the more entangled the electronic state. 

Our results support the idea that quantum information 
theory offers a powerful approach to the study of nonlin- 
ear complex system. At the transition point, the concur- 
rence is strongly affected, just as behaviors of the con- 



currence in the quantum phase transition point |22L \2c 
It would be more attractive to study entanglement be- 
haviors in other nonlinear complex systems such as two- 
dimensional quasicrystals and disorder systems. 
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